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EVALUATION OF AERODYNAMIC DERIVATIVES 


UNIVERSITY OF VIRGINIA WIND-TUNNEL COLD MAGNETIC 

BALANCE SYSTEM 

ABSTRACT 

The dynamic testing of a model in the University of 
Virginia cold magnetic balance wind-tunnel facility is expected 
to consist of measurements of the balance forces and moments, 
and the observation of the essentially six degree of freedom 
motion of the model. The aerodynamic derivatives of the model 
are to be evaluated from these observations. This paper is 
concerned with demonstrating the basic feasibility of extracting 
aerodynamic information from the observation of a model which 
is executing transient, complex, multi-degree of freedom motion. 
It is considered significant that, though the problem treated 
here involves only linear aerodynamics, the methods used are 
capable of handling a very large class of aerodynamic non- 
linearities. The basic considerations include the effect of 
noise in the data on the accuracy of the extracted information. 
For this purpose the motion of the projected first model, a 
cone, is calculated with "known" aerodynamic derivatives. This 
data is corrupted by various amounts of noise. The aerodynamic 
derivatives are then evaluated from the noisy motion by two 



methods: 1) the "Brute-force" method, and 2) the method 

of parametric differentiation. The "Brute-force" method 
treats the differential equations of motion as algebraic 
equations in the unknown aerodynamic derivatives and uses the 
method of least squares to average a large number of data 
points. In the method of parametric differentiation the 
equations of motion are considered in "aerodynamic derivative" 
space and the method of least squares is applied. Both methods 
extract the significant groupings of aerodynamic derivatives 
rather well and with only minor differences in accuracy. The 
relationships between noise level and the accuracy of the 
evaluated aerodynamic derivatives are presented. 
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A 


LIST OF SYMBOLS 


matrix of order n x n 
elements of matrix A 
matrix of order 1 x M 
elements of matrix B 

vector whose elements are the parameters 
3 C k /. 3(W/V t ) K = z , m 

3C k /3 (0.;d/V T ) 

3C k /3 <Wd/V 2 T ) 

3C K /3 (PdV/V 2 T ) 

3C r /3 (PVd 2 /V 3 T ) 

3D k / 3 (PRd 2 /V 2 T ) 

— (— ) 

V T 3t 

Base diameter 

Nondimensionalized moments of inertia, e.g. 
(A/(|pSd 3 ) ) 

Number of parameters 
Number of data points 
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P,Q,R 

p»q » r 
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Sj k ( T/C° ) 

SSR 
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u,v,w 


u, v,w 
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X,Y,Z 
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x,y ,z 


u 

p 
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LIST OP SYMBOLS (cont.) 

Angular velocities about body-fixed axes 

Pd/V T ,Qd/V T ,Rd/V T 

Characteristic area 

Parametric Coefficient 
Sum of square of error 
Time 

Translational velocities in body-fixed 
axes system 

Nondimensionalized velocities 

u/v T , v/v T , w/v T 

Tunnel speed 

Tunnel fixed orthogonal system of 
axes X positive downward 
k th element of X 

Nondimensional distance along (X,Y,Z) 
X/d, Y/d, Z/d 
Mass of model/ / ('|-P''S 1 d) 

Free stream air density 
Standard deviation 

Euler angles defining the orientation 
of body-fixed axes Cx, Cy, Cz (Cx being 
the axis of symmetry) in relation to 
tunnel-fixed axes 
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T Htv T /d) 

Subscripts and Superscripts 

i,j,k,l Integers 

o Initial conditions or values 
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Introduction 


For several years the development of magnetic balances to 
support models in wind tunnels has been in progress. The 
ability to support a model with no physical connections to it 
and the ability to arrange no model motion even at a lifting 
incidence, allow nearly free motion^or to influence and adjust 
the model motion to a significant extent is considered an 
important advance in wind tunnel techniques. Two different but 
complementary systems have been or are being developed. 

One of these, the MIT/French system, basically supports a 
cylindrical magnetic core and is effectively a five component 
balance. Roll control has been added but is essentially an 
independent subsystem. This system easily makes static tests 
of models in which the model is held at a fixed orientation and 
the forces and moments are measured. The system also can drive 
the model in simple oscillatory modes of motion to measure, for 
example, dynamic stability derivatives. Dynamic testing with 
the MIT/French system suffers from (1) the small magnitudes 
of the aerodynamic forces and moments caused by the motion in 
comparison with the forces and moments required to force the 
motion and (2) the complexities involved in dynamic calibrations 
of the magnetic balance system. 

In the other magnetic balance, the University of Virginia 



system, a magnetic sphere is supported and the system is 
basically a three component balance. By responding to trans- 
lation at frequencies below a somewhat adjustable "balance 
frequency", it holds the model in the test section. A model 
built around the magnetic sphere is nearly rotationally free 
at all frequencies and is essentially free in translation 
above the "balance frequency" (of the order of 10-15 hertz in 
the prototype) . This sort of operation is called Quasi-6 
Degree of Freedom operation and is considered to be the nearest 
thing to free-f light obtainable in a wind tunnel. The study 
of the dynamic stability of a model (a prime application of the 
U. Va. system) corresponds to measuring the balance forces and 
moments, observing the model motion, and inverting the equations 
of motion to obtain the aerodynamic forces and moments, static 
and dynamic. An obvious inconvenience with this system is that 
the motion of the model cannot be restricted to simple one or 
two degree of freedom oscillations. Thus the problem of ex- 
tracting the aerodynamic information from the experimental 
data involves the use of 6 degree of freedom motion data. The 
U. Va. prototype is designed specifically for the case of 
spinning axisymmetrical bodies (missile types) in supersonic 
flow. The extension to airplane types and subsonic flow is 
expected to be made in the not-too-distant future. 
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This paper describes and summarizes the results of a 
first attempt to demonstrate the feasibility of extracting 
aerodynamic information, specifically values of aerodynamic 
derivatives, from six degree of freedom motion data for a 
reasonably realistic approximation of a model in the prototype 
U. Va. system. Also of considerable interest is the relation- 
ship between noise in the data and accuracy of the extracted 
derivatives. The two methods of solution of the inverse 
problem ("Brute Force" {BF} and "Parametric Differentiation" 

{Par Dif f .})< which were used are sufficiently general that 
the inertial nonlinearities are retained and a quite large 
class of aerodynamic nonlinearities may be handled in principle 
as easily as linear aerodynamics. Only a transient motion 
case was considered. 

The general procedure was as follows. The projected first 
model for the U. Va. prototype system, a 15° cone, was chosen as 
the body and a set of "best" aerodynamic derivatives., was- 
established. With suitable approximations the full six degree 
of freedom equations of motion were numerically integrated to 
produce perfect motion data, which could be corrupted with 
various amounts of noise. The inertia parameters, balance 
forces and moments (if any) and the (noisy) data were used 
to invert the equations of motion to recover a set of aero- 
dynamic derivative. A comparison of the extracted derivatives 
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with those used to calculate the perfect data and the knowledge 
of the amount of noise added enables one to study the relation 
between noise level and accuracy of the extracted derivatives. 


Generation of Perfect Data 

For this study the base motion case is taken to be the 
projected first model, a 15° cone, in the U. Va. prototype 
cold magnetic balance wind tunnel system. The following 
assumptions and approximations are made: (1) Zero-gravity; 

(2) Zero-drag; (3) Zero rolling moments on the model; (4) Zero 
balance forces and moments on the model; (5) The mass distri- 
bution of the model is nearly axially symmetrical about the cone 

-2 

axis, but unbalanced so that i__ = i_ = i„ - 10 i_, {i- = i_; 

D £j r -D ID C 

i - i /6}; (6) The model is assumed to have a non-zero C , 

A o Z 


C,C ,C ,C,C,C , and C 
z z m m ' m_ ' m 

q Dw pv w q Dw pv 


(and, of course. 


the other corresponding lateral derivatives which follow from 
axial symmetry) [6, 7]; (7) The flight is at Mach 3; (8) The 

center of mass is on the symmetry axis at 0.6 of the length 
of the model from the nose. 

Subject to these assumptions, the full six-degree of 
freedom equations of motion are integrated by the Runge-Kutta 
4 step method. The translational equations are written in the 
tunnel fixed frame of reference, and the rotational equations 


9 



are written in the body fixed frame. Thus the perfect data 
consist;:' of sets of values of x, y, z, position of the center 
of mass in the tunnel fixed frame and ip, 0, <j> the conventional 
orientation angles. The initial conditions are all zero except 
i(j(0) = 0.1, 0(0) = 0.15, and ^4 (0) = 0.016 5 a roll rate close 
to roll pitch resonance for this model. The motion of a model 
in the U. Va. facility is expected to__be observed as a continuous 
record in time, which may then^ -reduced -to di s ere ^ data 
points as dense as desired. About 4 cycles in 0 and ip , 100 
points, are used for the positional (x,y,z ,i|>,6,<f») perfect data. 

Inverse Problem 

The inverse problem is solved by two methods: The "Brute- 

force" method and the method of parametric differentiation. 

The aerodynamic axial force and rolling moment are assumed to be 
zero in the generation of the perfect data and are not considered 
in the inverse problem. The aerodynamic derivatives are evalu- 
ated in the body-fixed coordinate system. 

(1) The "Brute-force" (BF) Method 

The BF method of solution uses the discrete positional 
data, numerically calculates velocities and accelerations and 
inserts them, along with other known quantities, into the 
equations of motion. The equations of motion are thus treated 
as linear algebraic equations in the 
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der ivatives .The full nonlinear character of the inertial terms 
in the equations of motion is retained in the solution. 

The model being axisymmetric , only the z force (in body 
frame) and pitching moment equations of motion are considered 
in the inverse problem. Furthermore, the force and moment 
equations are independent of one another. 

The pitching moment equation, for example, is written as 
an algebraic equation in the unknown aerodynamic derivatives 
as follows: 


m 


w 


w.+C q.+C Dw.+C (pv) . = i„ Dq . 
1 m_ m Dw 1 m 1 B 


pv 


( 1 ) 


- (i c - i A ) (rp) ± + i E (p 2 ~ r 2 ) ± - i p (qr + Dp) ± + i D (qp -Dr) ± 
i = 1, 2, , N, 


where N is the number of data points. This can be written as 
follows : 

n 

j=l a ij C j = b i i = 1, , N (2) 


where c.'s are the moment derivatives, 
of motion, and b^'s the inertial terms 
an error equation 


a. .'s are the kinematics 
Rewriting Eq. (2) as 


N n 

SSR = E (E a. . c . - b.) 2 (3) 

±=I . j=l ^ 3 
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and minimizing SSR with respect to c^, the aerodynamic deri- 
vatives are obtained. The z force equation is treated in the 
same fashion. 

The kinematics of motion, u, v, w, p, q, r, Du,,. Dv, Dw, 

Dp, Dq, Dr are obtained from the positional data in the 
following manner: In case of noisy data, the positional 

data at every instant are smoothed by using a quintic power 
series and 9-equidistant points. The velocities and accelera- 
tions in the tunnel fixed reference frame are obtained by 
numerically differentiating 5-equidistant points of the 
smoothed translational positional data. The Euler angular 
velocities and accelerations are obtained in a similar fashion. 
These are transformed to the body axis reference frame to 
obtain the necessary kinematics of motion. 

(2) The Method of Parametric Differentiation 

The motion of the model is considered as being in 
a parametric space of aerodynamic derivatives and the initial 
conditions of motion. The six equations of motion are used 
and are written as a set of first order nonlinear differential 
equations of the form 


DX = F (X, C, t) (4) 

where X is the set of twelve (12) positions and velocities and 
C is the twelve initial conditions, X(o), plus the aerodynamic 
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derivatives to be evaluated. For a given set: of c° the 
solution is of the form 


X = X (t, C ) 


The change in X due to small changes in the parameters C is 
given by the sensitivity coefficients, defined thus: 


S k (x, C°) = lim 
^ Ac^-0 


k . o o o , • o, k . _,o. 

x (x , c lf c 2 , — c. .+. .Ac ..,. - — , c^) - x (t, C ) 


3x k (x , C°) 


Differentiating Eq. (4) with respect to c^ it follows that 

DS k (t, C°) =' (x, C°) + (7) 

3 - 3x £ 3 - 

These are the parametric equations, which are a set of 
linear differential equations with variable coefficients. 

The solution of Eq. (4) is expressed in terms of Eq. (5) 
and Eq. (6) in the form of a truncated Taylor series: 

x k (x, C) = x k (t, C°) + Qe 11 S k ( , C°) Ac... (8) 

i=l 1 1 

In the present case the parameters C are ordered sequentially 
in the following way: 


x , y , z , ip , 0 , <f> , Dx , Dy , Dz , Dip , D6 , D<j> , 
o J o o o o o o 2 o o o o o 


C , C , C , C , C , C , C ? C 
m m m m„ z z z., * 

W q pv Dw w q pv,. Dw 
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With this sequence it follows from consideration of Eq. (8) 
that 


Sj k (0, C°) = 1 if j = k 
= 0 if j ^ k 


(9) 


Equations (7) and (9) constitute the basis for the para- 
metric study of the motion governed by Eq. (4) . The observed, 
positional data, in tabulated form, can be denoted by 

X,. 

> k = 1 , — 6 and i = 1 , N 

k o 

Let x^ (t, C ) be the values of the calculated motion obtained 
by solving Eq. (4) with estimated values of C°. The sum of 
squares of the error between the observed and calculated 
positions is given by the following: 


N 6 , 

SSR = Z ’Z [ e x. ' 
i=l k=l 1 


(x i N (t, C°) + Z S k . Ac .)] 2 

1 j=l J1 3 


Minimizing SSR with respect to Ac^ one gets the following: 


where 


A AC = B 


N 6 k k 
a. = Z Z S K . S,. 
3* i=l k=l 31 *1 


N 6 , 

b . = Z^ z {x . : <- e x . ) S k 


i=l k=l 


jk 
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This yields the correction to the assumed C° as 

AC = A -1 B 

The new set of C° risj used in Eq. (4) and the procedure con- 
tinued until convergence in the SSR is obtained. The" final set C is 
defined to be the "true" C. 

The calculated motion depends on the initial estimation 
of the aerodynamic derivatives and the initial conditions 
of the motion. The latter are estimated by obtaining the 
smooth positional data and the velocities by numerical differenti- 
ation of the first few data points by the methods used in the 
"Brute-force" method. 

Modelling Problem 

The modelling problem is that of determining a "best" set 
of aerodynamic derivatives for a body from its observed motion. 

It is conceivable that an important derivative could be missed, 
if the aerodynamic characteristics of the body are sufficiently 
strange and unknown. 

A very preliminary attack on the general modelling problem, 
as well as a basic verification of the two chosen inversion 
methods, is accomplished in the following way. One assumes the 
perfect data are real, and of course very accurate, motion data 
for a body, the aerodynamic derivatives. of which are completely 
unknown. The problem then is to see if one can construct a 
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logical argument to determine a "best" set of aerodynamic 
derivatives for that data subject to the ignorance assumption. 
(Of course, in the present case, one may also make comparisons 
and draw conclusions using the "known" derivatives for that 
data . ) 

The process of interpretation of the results is affected 
by the following considerations. The aerodynamic pitching 
moment for the chosen model is 


C = C w + C q + C Dw + C pv, 
mm m ^ m_ m 

w q Dw pv 


By adding and subtracting suitable terms the same pitching 
moment may be written in three alternate forms: 


C = C w + (C + C ) q + (C - C )pv + C (Dw - q + pv) 

m m m itl ^ m „ nu ^ 

w q Dw pv Dw Dw 

= C w + (Cr + C ) q + (C - C ) Dw + C (Dw - q + pv) 

m mm m m 

w q pv Dw pv pv 


= C w + (C + C ) pv + (C + C - )Dw - C (Dw - q + pv) 
m mm m_ m m 

w q pv Dw q q 


For a low lift configuration, such as a 15° cone, one expects 
the combination (Dw-q+pv) to be small. If Dw - q + pv is a 
one order smaller quantity than q, Dw, or pv, then one may expect 
only two damping constants to be well determined. The second 
and third terms on the right sides are three alternate forms of 
these two groupings. Thus one can expect that the values of 
these damping derivative groupings'' to; be of more significance than 
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the individual derivatives. Linear combinations of any pair 
will give any other pair, hence there is little significance 
in which pair of groupings one uses. Exactly similar consid- 
erations hold for the significant groupings of force damping 
derivatives . 


(1) The "Brute-force" Method 

The number of aerodynamic derivatives, n, chosen in 
Eq. (3) is arbitrary. This is true also for the corresponding 
force equation. Table 1 shows the result of choosing n = 2, 

3, 4, 5 and extracting 2, 3, 4, 5 individual force and moment 
derivatives. For n = 6, adding C and C to the set of 

if w;; 

derivatives, the matrix A was singular and no solution was 

found. Table I also shows the important groupings which follow 

from the fact that the quantity Dw + pv - q is very small. 

An inspection of Table I from the "modelling Problem" 

point of view (derivatives unknown) is quite informative. For 

these four sets of extracted derivatives, one observes that the 

values of C and C are nearly constant. C , C , C , C 

m z z in z _ m_ 

w w q q Dw Dw 

seem to oscillate with increasing n. The two groupings of the 

damping derivatives remain nearly constant for n 3 . It is 

apparent that, when two or more damping derivatives are extracted, 

(C + C ) and (C - C ) , as well as the corresponding 
m m m_ m 

q pv Dw pv 

force groupings are nearly constant. 

It seems clear that this particular body has two important 


17 



moment damping derivatives and two important force damping 
derivatives which show in the groupings regardless of how many 
derivatives are extracted. The fact that the force derivatives 


behave nearly as well as the moment derivatives is somewhat 
surprising . 

An inspection of Table I including the actual values and 
the percent errors verifies the above "modelling problem" 
conclusions. Large errors in the individual derivatives can 
occur, e.g., n = 4, with excellent values for the two inde- 
pendent groupings. One can conclude that for this body in 
this motion state that when two or more force and moment damping 
derivatives are included in the set successfully extracted from 

accurate data, the groups (C + C ) , (C - C ) and the 

q pv Dw pv 

corresponding groupings of force derivatives are given with 
good accuracy. 

(2) The Method of Parametric Differentiation 
The same perfect data are sused for extraction of derivatives 
by the method of parametric differentiation. Increasing the 
number of parameters from 15 to 20 corresponds to taking 


and in order C 


C in 
z Dw 


the inverse problem. The ordering of the aerodynamic derivatives 


is arbitrary but in the work reported here that order was fixed. 


For all the cases in Table II, except that in the last column. 


the starting values were the actual values. For the last case 
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Inversion of the Perfect Data by the Brute Force (BF) Method 





INVERSION OF THE PERFECT DATA BY THE METHOD 





0.4 


(n = 19, 5 iterations) the initial estimates were C = 

m 
w 

and C = -1.8 and the rest were zero. For the missing n = 16 
z 

w 

case (four- moment derivatives and no force derivatives) the 
program gave numbers which were unrealistic (magnitudes too 

3 

large by 10 or more) . 

For this method an obvious and natural criterion for the 
"best" set of aerodynamic derivatives is: Minimum SSR corresponds 

to the "best" set of derivatives. It is probably no real dis- 
advantage that this criterion loses its significance when SSR 
becomes sufficiently small. The Modelling Problem point of 
view, i.e., disregarding actual values and percent errors, again 
may be adopted in viewing the results presented in Table II. 

A first conclusion is that a few iterations are very 
helpful in reducing SSR as shown by the n = 15 and n = 19 cases. 
Presumably a few iterations for n = 17 and 18 would have reduced 
SSR significantly. Secondly, the moment damping groupings 
change much less as n is changed than the individual moment 
damping derivatives. A similar conclusion about the force 
damping derivatives cannot be made since the force derivatives 
occur in so few cases. 

The examination of the entire sequence suggests: (1) the 

translational motion is important in this motion as evidenced 
by n = 16 not working and the general reduction in SSR as force 
derivatives are added, (2) there are two important moment damping 
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constants corresponding to the groupings, and (3) the sharp 
reduction in SSR in going from n = 18 to n = 19 suggests that 
there are (perhaps only) two important force damping constants. 

Again, adding the actual values of the derivatives and 
the percent errors gives quantitative evidence to the above 
conclusions. The accuracy of the two groupings of force 
damping derivatives in case n = 19 is confirmed. The inclusion 
of translational motion in the inversion problem does increase 
the accuracy with which the moment damping groupings are 
extracted. 

There is one interesting item in the comparison of the two 
methods. There is only one common case tried by each. Brute 
Force, n = 4 (Table I) . corresponds to extracting the same 
eight non-zero derivatives used to generate the perfect data. 

It was successful in that values of the important force and 
moment groupings came out with fair accuracy, though all the 
percent errors were considerably worse than for the adjacent 
cases (n = 3,5), The Parametric Differentiation case n = 20 
corresponds to extracting the same set of derivatives, and it 
did not work. It is tempting to conclude that Brute Force is 
a somewhat more stable numerical process than parametric 
differentiation. 

Influence of Noise 

In order to investigate the influence of noise in the data 
on the accuracy with which aerodynamic derivatives may be 




extracted, one case of each of the two inversion methods was 
chosen. These were Brute Force n = 3 and Parametric Differentia- 
tion n = 19, which seemed to give best results with the perfect 
data . 

Noisy data .were simulated by adding to the data a set of 
pseudo-random numbers with zero mean and a chosen standard 
deviation. The noise levels in the translational positions 
x,y, and z were chosen to be three times larger than those in 
the angular positions ip and 0. This corresponds crudely to 
some notion of how accurately translational and rotational 
positions may be measured. Noise was not added to the roll 
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angle <j>*. At each of the noise levels five (5) experiments 
(independent sets of noise) were done. Different or inde- 
pendent sets of pseudo-random numbers with zero mean and the 
same standard deviation were generated. In the case of "no 
noise" the five experiments correspond to 5 sets of slightly 
altered initial conditions used to calculate the perfect data. 
The results are shown in Tables III and IV. ParDiff was 
done with 3 iterations. 


♦Originally, the noise added to the perfect data was of the 
form 


X = X (1+Y ) 

-noisy -perfect -noise 

i.e., relative noise. For a reason which was not clear at 
the time, the error in the extracted aerodynamic derivatives 
was especially sensitive to (relative) noise in <t>. The change 
from relative to absolute noise 

X = X + Y 

-noisy -perfect -noise 

was made because the latter seems to correspond better to 
reality. After approximately one half the calculations were 
done a case with absolute noise added to <J> showed that the error 
was quite insensitive to (absolute) noise in cf>. For con- 
sistency (and the lack of computing funds) the calculations were 

completed with no noise in <f> . It is now clear that the secular 

* 

character of <j>, <j) is nearly constant, causes a small relative 
noise to correspond to a quite large effective absolute noise. 
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PARAMETRIC DIFFERENTIATION, n = 19 (PAR DIFF) INVERSION OF NOISY DATA 



Tables III and IV present, for. each derivative and 
important derivative grouping at each noise level, the mean of 
the sample of five numerical experiments and the standard de- 
viation for the sample. In almost all cases the standard 
deviation is considerably larger than the error in the mean. 
Further, the errors in the mean show no definite trend with 
noise level, while the standard deviation shows a definite 
increase as the noise level increases. Therefore, the standard 
deviation, rather than error in the mean, is related to probable 
error in the values of extracted derivatives. It is to be noted 
that the confidence level is low because the number of experi- 
ments in a sample is small. The largest noise level, cw = 

- ° z = -009 correspond to o being about 10% of the amplitude 
of a typical position variable. 

It is not surprising to note, that the static derivatives 
are extracted with greater accuracy than are the dynamic deri- 
vatives. It may be surprising to note that Parametric Differenti- 
ation recovered C out of the noisy data with about twice the 

z 

w 

accuracy as it recovered C . Table V presents the standard 

w 

deviation in percent of the static derivatives as a function of 
the noise level for the two methods of inversion. 

Figure 1 shows how the standard deviation in percent 
of the force and moment damping derivative groupings varies with 
the noise level for both methods. An overall view of Figure 1 


27 



• 

• 

• 

• 

Q 

a 




o 

o 

O 

O 

X 

o 




o 

o 

o 

o 

II 

M 




VD 

CO 

h- 1 

o 


c o 






Cn 

'•J 

Q 

w 







Cn 

v; 









n 

t- 1 








Q 

pi 








N 

< 









M 









f 




oo 

H-* 

!-* 







• 

• 

• 

• 

o 

dd 





Co 

O 

NO 

N 

pd 




c> 

Cn 

Co 

Cn 

s, 

a 




CO 

4S 

Ov 

4N 


H 









w 


CO 









H 







pj 









o 









pd 


0 







n 


> 







w 


g 


NO 

h-» 



o 


Hi 

a 


• 

• 

• 

• 

3 


O 

© 


CO 

M 

Cn 

h-» 

« 


H 

w 


o 

h- 1 

Ov 

CN 




< 


oo 

OO 

O 




CO 

H 








H 

> 








> 

H 








H 

M 

H 







H 

O 

> 







n 

a 

dd 









a 







© 

M 

M 







w 

a 








pj 


< 







H 









< 









> 

< 








H 

cn 


f— 1 




n 

a 

M 



• 

• 

• 

• 

N 

> 

< 

a 


o 

NO 

O 

O 



M 

o 


Cn 

NO 

VO 

VO 


3> 

CO 

H 



VO 

Ov 

Ov 

o 

2* 


CO 






M 

pj 


PI 






PI 

H 








>n 



a 






pi 

M 


pj 






2 

O 


<3 






w 



W 






a 



a 






H 









M 









> 









H 





NO 

M 



n h 





• 

• 

• 

• 

3 O 





N3 

h-» 


CO 

s a 





M 

CO 

4N 

Cn 






-O 

Cn 

NO 

O 









28 



shows that Par Diff gives better results. More specifically. 

Par Diff recovers the force groupings about four times better 

than BF. Par Diff extracts (C + C ) and (C - C ) 

m m n iru m 

q Dw Dw pv 

equally well, while BF extracts the former about three times 
better than the latter and the former not quite as well as 
does Par. Diff. 

The computing requirements for the two methods are signi- 
ficantly different. Doing 3 iterations. Par Diff requires about 
25 times as much computing time as BF . 

Conclusions 

(1) There is no particular difficulty in extracting aero- 
dynamic information of reasonable accuracy from reasonably 
noisy data describing the complicated, Quasi-Six degree of 
freedom motion of a 15° cone in a Mach 3 tunnel. This result 
is demonstrated by a quite simple and straightforward inversion 
technique (BF) and by a quite sophisticated and complex inversion 
technique (Par Diff) . The interesting implication is that this 
conclusion is not unique to this particular motion case. 

(2) Both of the inversion methods used are inherently 
general in that they do not depend upon linearization. BF is a 
linear algebraic problem as long as the aerodynamic forces and 
moments may be written in terms of the derivative coefficients. 
Par Diff involves a numerical integration of the equations of 
motion; the inertial nonlinearities present no difficulties, but 
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it is conceivable that enough higher order terras in the 
expressions for the aerodynamic forces and moments , could be im- 
portant and could cause excessive integration difficulties. 
Certainly most important aerodynamic non-linearities could 
be handled by Par Diff. The interesting possibility is that 
these methods, and perhaps others, coupled with magnetic 
balance wind tunnel systems provide a practical way of ex- 
perimentally studying nonlinear aerodynamics of suitable 
models . 

(3) The work reported here is only a first step in 
the attempt to understand the aerodynamic information extraction 
problem for the complicatedn;. guasi-s ix degree of freedom case. 
The obvious extensions which ought to be done involve: various 

noise levels in individual observables and their influence in 
combination, different models and flight conditions, and active 
balance forces and moments. 
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Figure I. Standard deviation in % of the extracted derivative groupings as 

a function of noise level. The highest noise corresponds to a* being 
about 10% of the amplitude of the typical variable. 


